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■ PACS. 47.65+a - Magnetohydrodynamic. 

PACS. 47.27.GS- Isotropic Turbulence. 
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We investigate the time evolution of two different (GOY-like) shell models which have been 
recently proposed to describe the gross features of MHD turbulence. We see that, even if they 
are formally of the same type sharing with MHD equations quadratic couplings and conserved 
quantities, differences exist which are related to conserved quantities. 
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One of the striking features of turbulence in fluid flows is the existence of an energy cascade 
from large to small scales, which is usually described in the framework of the Richardson's 
picture [[l]. Some efforts to describe the cascade through toy models or phenomenological 
considerations have been made in the last years 0, aimed to understand chaotic behaviour, 
energy spectra, intermittency, etc. Among others, Shell Models can be built up by dividing the 
wavevectors space k into logarithmically spaced shells A™ < k/ko < X n+1 (A is usually taken 
equal to 2), each shell being characterized by a discrete wavevector k n — X n ko (n = 0, 1, N), 
and by one single dynamical variable whose evolution is representative of the dynamics of fields 
belonging to the shell. The equations for these variables can be written by assuming quadratic 
nonlinear couplings between shells (as the original equations), and some of the coupling 
coefficients can be fixed by invoking the conservation of invariants in the dissipationless and 
force-free case. Up to now the so called GOY-model ^ ^ has been considered as the most 
representative to reproduce the turbulent features in fluid flows. 

Shell Models for MHD turbulence have been introduced in the last years pf , and in particular 
the properties of GOY-like models have been investigated by Biskamp |6| and more recently 
by Frick and Sokoloff (7) (herein after referred to as B model and FS model respectively). 
Both models can be obtained as particular cases of the following set of equations (v n and b n 
represent respectively the velocity and the magnetic field in dimensionless units) 

= -vk 2 n u n - v'k~ 2 u n + ik n ^[u n+1 u n+2 - b n+ ib n+2 ) 

- &n-l&n+i) j2~{u n -2U n -i - b n - 2 b n -l)} + fn (1) 
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= -r\k 2 n b n + ife n |(l - 5 - <5 m )(zt n+ i& n+2 - b n+ iu n+2 ) 
6 1 — 5 i * 

H — ^-(w n _l6„+l — 6 n _i« n+ i) H ^ ("n-J^n-l ~ &n-2Wn-l)j + ffn (2) 

or, in terms of the complex Elsasser variables Z„(t) = v n (t) ± & n (£), particularly useful in 
some solar-wind applications, 



dZ- 



where 



(// -v + klZ±-v-k 2 n Z*--k- 2 Z+--k- 2 Z-+ik n T±* + ft (3) 



± n — \ 2 Zj n+l Z/ ?i+2 i 2 Z/ n+l z 'n+2 



S rn - 5 ± T _ 5 + S m T ± 

2^2 ^n-l Z ri-2 ^ Vl Vj f l4) 
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Here v ± = (u ± r?)/2, being the kinematic viscosity and 77 the resistivity, 
eq. (jl]), is a drag term specific to 2D cases (see below), / ? f = (/„ ± g n )/2 are external 
driving forces, 6 and <5 m are real coupling coefficients to be determined. In the inviscid 
unforced limit, equations (Rh conserve both pseudoenergies E ± {t) = (1/4) \Z± (t)\ for 
any value of 8 and S m (the sum is extended to all the shells), which corresponds to the 
conservation of both the total energy E = E + + E~ = (1/2) ^2 n (\v n (t)\ + \b n (t)\ ) and the 
cross-helicity he = E + — E~ = Re{v n b* n ). The values of 5 and 5 m are fixed by imposing 
the conservation of another quantity which is the magnetic helicity in 3D or, in 2D, the mean 
square magnetic potential [^j. In analogy with the fluid case Q we can define a generalized 
quantity as H^\t) = X)n( s ig n (^ ~ ^)) n \bn(t)\ 2 /k% whose conservation implies 5=1 — A~ Q , 
§ m = \~ a /(l + \- a ) for S < 1, < 6 rn < 1 and 6 = l + A - ", 5 m = -\- a /{l - X a ) 
for 5 > 1, S m < 0, 6 m > 1. Thus two classes of MHD GOY models can be defined with 
respect to the values of 5: 3D-like models for S < 1, where is not positive definite; 

2D-like models where 5 > 1 and Hg is positive definite. This situation strongly resembles 
what happens in the hydrodynamic case where 2D-like (S > 1) and 3D-like (5 < 1) models 
are conventionally distinguished with respect to a second generalized conserved quantity 
H< k\t) = ]T n (sign(<5 - l)) n K \v n {t)\ 2 . Here the 3D and 2D cases are recovered for a — 1, 2 
where the ideal invariants are identified respectively with kinetic helicity and enstrophy. It 
should be noted that, although the hydrodynamic invariants are not conserved in the magnetic 
case, the equations which link a and 5 are exactly the same for hydrodynamic and MHD 
models. Thus, once fixed a and 5, it is a simple matter to find out which GOY model the 
MHD GOY one reduces to when b n = 0. To summarize we have that (with A = 2) the FS 
model for the 3D case is recovered for a = 1, 5 = 1/2, S m = 1/3 and reduces to the usual 
3D GOY model for b n = 0. The B 3D model is actually a 2D-like model. It is obtained for 
a = 1, 5 = 3/2, S m = — 1 and reduces to a 2D-like GOY model that conserves a quantity 
which has the same dimensions as kinetic helicity but is positive definite. The 2D FS and B 
models coincide, they are recovered for a = 2, 5 = 5/4, S m = —1/3 and reduce to the usual 
2D GOY model for b n = 0. 

Let us now briefly review the main results found in The authors investigate the problem 
of the magnetic field generation in a free-decaying turbulence, thus showing that: 1) in the 
31? case magnetic energy grows and reaches a value comparable with the kinetic one, in a 
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Fig. 1. - FS model: log 10 (|v„| 2 ) (diamonds) and log 10 (|6„| 2 ) (lines) versus log 10 fc n . The 
averages of \b n \ 2 are made over intervals of 3 large scale turnover times. Time proceeds 
upwards. The kinetic spectrum is averaged over 30 large scale turnover times. The straight 
line has slope -2/3. N=24, v = rj = 1CT 8 , u' = 0. 

Fig. 2. - B model: log 10 ( \v n \ 2 ) (diamonds) and log 10 ( \b n \ 2 ) (lines) versus log 10 k n . Averages 
are made over intervals of 100 large scale turnover times. Time proceeds downwards. The 
kinetic spectrum is only shown for the last interval. The straight line has slope —4/3, see text 
for explanation. N=33, v = lO" 16 , r\ = 0.5 • 10" 9 , v' = 1. 



way that the magnetic field growth is unbounded in the kinematic case; 2) in the 2D case 
magnetic energy slowly decays in the nonlinear as well as in the kinematic case. These results 
have been interpreted as a 3D "turbulent dynamo effect" Q| and seem to be in agreement with 
well-known results by which dynamo effect is not possible in two dimensions |To|j . Morover 
the authors investigate the spectral properties of the model in the stationary forced case, and 
they find that Kolmogorov spectral properties are only established when the cross-correlation 
is smaller than the energy of the system. In this paper we are going to compare the results 
obtained from the 3D B model and the 3D FS model as far as the "dynamo-like effect" is 
concerned. We started from a well developed turbulent velocity field and injected a seed of 
magnetic field looking at the growth of the magnetic spectra. System is forced on the shell 
n = 4 (fco = 1), setting = = (l + «) 10 -3 , which corresponds to only inject kinetic energy 
at large scales. Method of integration was a modified fourth order Runge-Kutta scheme, with 
a time step 10 -4 . In fig. 1 we report log 10 ( |& n | 2 ) and log 10 ( \v n \ 2 ) versus log 10 k n for the FS 
model. Angular brackets ( ) stand for time averages. It can be seen that the magnetic energy 
grows rapidly in time and forms a spectrum where the amplitude of the various modes is, at 
small scale, of the same order as the kinetic energy spectrum. Besides the spectral index is 
very close to /c~ 2 / 3 . For a comparison we integrated the 3D B model and it can be seen (fig. 
2) that a magnetic spectrum is formed, but it slowly decays in time. Notice that, because of 
the smallness of b n , its back-reaction on the velocity field is negligible, thus the kinematic part 
of the model evolves independently from the magnetic one. Then the scaling \v n \ 2 ~ k n 4 ^ 3 
follows, in agreement with || where a cascade of generalized-enstrophy is predicted for 2D-like 
hydrodynamic GOY models when a < 2. The question now arises whether it is correct the 
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Fig. 3. - Ideal case: kinetic energy (continuous line) and magnetic energy (dashed line) versus 
time for the FS model; magnetic energy (dot-dashed line) versus time for the B model. 



interpretation of the growth of the magnetic field in the FS model as the corresponding dynamo 
effect expected in the real 3D MHD. First of all it should be noted that in the kinematic case 
an analogy with the vorticity equation predicts the following relations between velocity and 
magnetic energy spectra ||: \v n \ 2 ~ fc~ a , \b n \ 2 ~ k 2 ~ a , so that if a = 2/3 it follows a 
magnetic energy spectrum growing with k. The kinematic case corresponds to the first stage 
of growth of our simulation where this behaviour is sometimes visible, at least qualitatively. 
Note however that the averages are made on very small time intervals because of the rapid 
growth of the magnetic energy. A similar, much more pronounced behaviour is found for the B 
model as well. As regards the final state of evolution for the FS model, it should be reminded 
that equipartition of magnetic and kinetic spectra at small scales is generally attributed to the 
Alfven effect, that is the propagation in opposite directions of turbulent eddies Z+, Z~ along 
the largest scale magnetic field (see for example ||). Note however that, as already remarked 
in Q for another kind of shell models, equipartition is reached even if the Alfvenic terms are 
not included. The reader is referred to Q for a thorough discussion concerning the inclusion 
of Alfvenic terms in shell models. 

Actually we want to stress that: 1) The sign of the third ideal invariant seems to play a 
crucial role as far as the growth of small magnetic fields is concerned; 2) The stationary final 
state of evolution for the FS model does not represent a Kolmogorov-like turbulence. In fact 
let us consider the ideal evolution of the model dZ^ jdt = ik n T^*- We can build up the phase 
space S of dimension D = AN , by using the Elsasser variables as axis, so that a point in S 
represents the system at a given time. A careful analysis of (||) shows that there exist some 
subspaces I c S of dimension D = 2N which remain invariant under the time evolution [ prf . 
More formally, let y(0) — (v n ,b n ) be a set of initial conditions such that y(0) £ I, I is time 
invariant if the flow T*, representing the time evolution operator in S, leaves / invariant, that 
is T'[y(0)] = y(t) £ I. The kinetic subspace K C S, defined by y(0) = (v n ,0) is obviously 
the usual fluid GOY model. Further subspaces are the Alfvenic subspaces defined by 
y(0) = (v n ,±v n ), say Z£ ^ and Z~ = (or vice versa). Each initial condition in these 
subspaces is actually a fixed point of the system. We studied the properties of stability of 
K and A^. Following plj ], let us define for each / the orthogonal complement P, namely 
S = I® P. Let us then decompose the solution as y{t) = (yi n t(t), y e xt(t)) where the subscripts 
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refer to the / and P subspaces respectively. Finally we can define the energies Ei nt = \\yi n t\\ 2 

and E ext = \\y ex t\\ 2 - Note that the distance of a point y = {yi n t,y e xt) from the subspace / 

is d — min \\y — y\\ = \\y e xt\\- Then E ext represents the square of the distance of the solution 
yei 

from the invariant subspace. At time t = 0, E ext = eEi nt (e <C 1) represents the energy of the 
perturbation. Since the total energy is constant in the ideal case, two extreme situations can 
arise: 1) The external energy remains of the same order of its initial value, that is the solution 
is trapped near / which is then a stable subspace; 2) The external energy assumes values of 
the same order as the internal energy, that is the solution is repelled away from the subspace 
which is then unstable. Since the external and internal energies for the Alfvenic subspaces are 
nothing but the pseudoenergies E + and E~ , which are ideal invariants, the Alfvenic subspaces 
are stable. As regards the kinetic subspace, Ei nt and E ext represent respectively the kinetic 
and magnetic energies. Looking at the numerical solutions of the ideal model (fig. 3) we can 
see the difference in the stability properties between the B model and the FS model. In the 
first case the external energy remains approximately constant, while in the second case the 
system fills up immediately all the available phase space. This striking difference is entirely 
due to the nonlinear term, and in fact must be ascribed to the differences in sign of the third 
invariant. The effect of the unstable subspace, which pushes away the solutions, is what in 
ref. @ is called "turbulent dynamo effect". When we introduce constant kinetic forcing and 
viscosity, the stable subspaces become attractors, so that the B model is attracted by its own 
kinetic subspace, while the FS model is attracted towards one of the Alfvenic subspaces. The 
dynamics of the FS model is strongly dominated by this behavior. Long runs show that the 
system evolves inevitably towards a "dynamical alignment" in which v n — ±6 ra as can be seen 
in fig. 4 where we plot Re{vj) and Re{bi) versus time. It is clearly visible how strongly velocity 
and magnetic field are correlated. In fig. || we plot, for the FS model, the time evolution of 
the reduced cross-helicity (E + — E~)/ (E + + E~), which is a global measure of the dynamical 
alignment. It can be seen that even from an initial value he = it grows towards unity. 
Due to the particular form of the nonlinear interactions in MHD ml), the nonlinear transfer 
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towards the smaller scales tends to be stopped and v n and b n spectra become much steeper 
than that reported in fig. 1. In this case the turbulent statistical properties are not so clear as 
in a Kolmogorov-like turbulence. The observed behaviour of the FS model might be avoided 
with a careful choice of the external driving forces. A work in this perspective is actually in 
progress. 

*** 
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who gave us the preprint of his paper and to the referees whose comments have improved the 
final version of the paper. 



REFERENCES 

[1] Frisch U., Turbulence: the legacy of A. N. Kolmogorov, (Cambridge University Press) 1995. 

[2] Bohr T., Jensen M.H., Paladin G. and Vulpiani A., Dynamical systems approach to turbu- 
lence, (Cambridge University Press) 1998. 

[3] Gledzer E.B., Sov. Phys. Dokl., 18 (1973) 216; YAMADA M. and Ohkitani K., J. Phys. Soc. 
Jap., 56 (1987) 4210. 

[4] Jensen M.H., Paladin G. and Vulpiani A., Phys. Rev. A, 43 (1991) 798; Biferale L., 

Lambert A., Lima R. and Paladin G., Physica D, 80 (1995) 105. 
[5] Gloaguen C, Leorat J., A. Pouquet and Grappin R., Physica D, 17 (1985) 154; Grappin 

R., Leorat J. and Pouquet A., J. Phys. (Paris), 47 (1986) 1127; Carbone V., Europhys. 

Lett., 27 (1994) 581; Frick P., Magnetohydrodynamic, 20 (1984) 262; Carbone V. and Veltri 

P., Geophys. Astrophys. Fluid Dyn., 52 (1990) 153. 
[6] Biskamp D., Phys. Rev. E, 50 (1994) 2702; BlSKAMP D., Chaos, Solitons & Fractals, 5 (1995) 

1779. 

[7] Frick P. and Sokoloff D., Phys. Rev. E, 57 (1998) 4155. 

[8] Biskamp D., Nonlinear Magnetohydrodynamics, (Cambridge University Press) 1993. 

[9] P. D. Ditlevsen and I. A. Mogensen, Phys. Rev. E, 53 (1996) 4785. 
[10] Zeldovich Ya.B., Sov. Phys. JETP, 4 (1957) 460. 
[11] Carbone V. and Veltri P., Astron. Astrophys., 259 (1992) 359. 



